Elementary vortex pinning potential in a chiral p-wave superconductor 
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The elementary vortex pinning potential is studied in a 
chiral p-wave superconductor with a pairing d = z(k x ± ik y ) 
on the basis of the quasiclassical theory of superconductiv- 
ity. An analytical investigation and numerical results are pre- 
sented to show that the vortex pinning potential is dependent 
on whether the vorticity and chirality are parallel or antipar- 
allel. Mutual cancellation of the vorticity and chirality around 
a vortex is physically crucial to the effect of the pinning center 
inside the vortex core. 

PACS numbers: 74.60.Ge, 74.60.Ec, 74.70.Pq, 74.70.Tx 



Much attention has been focused on the vortex pinning 
in type-II superconductors. The vortex pinning plays an 
important role on various vortex-related quantities and 
phenomena such as the critical current and the hysteresis 
of the magnetization in superconductors under magnetic 
fields. The characteristics of the vortex-related phenom- 
ena are of particular interest in unconventional supercon- 
ductors with multiple components of the superconduct- 
ing order parameter £1 In such superconductors, multiple 
states of superconducting order can coexist. Accordingly 
there appear multiple kinds of vortex structure, where 
the nature of the vortex pinning can be dependent on 
the microscopies of the superconducting order. 

One of the superconductivity with multiple compo- 
nents of the order parameter is the chiral p-wave one, 
d = z(k x ± iky), which is composed of two degenerate 
pairing states k x and k y and breaks the time-reversal 
symmetry. This superfluid 3 He-A type of chiral p-wave 
pairing state has been anticipated in a layered ruthenate 
superconductor S^RuO^d While the identification of the 
genuine superconducting parang of this material is still 
open to further discussionjj a that chiral p-wave pairing 
has the simplest and essential form and has attracted a 
great deal of attention. The ^oxtices for Sr2Ru04 have 
been investigated intensivelyLTtEl In the context of the 
vortex pinning, we will see a rich physics contained in 
that chiral p-wave pairing state. 

In this paper, we investigate the elementary vortex pin- 
ning potential in the chiral p-wavc superconductor with 
the pairing d = z(fc x ± ik v ). A point-like pinning cen- 
ter and a single vortex with vorticity perpendicular to 
conduction layers in a layered superconductor are consid- 
ered. We show that the vortex pinning potential depends 
on the sense of the chirality of the Cooper pairs relative 
to the vorticity of the vortex in the chiral p-wave su- 
perconductor. First we analytically discuss the interplay 
between the chirality and vorticity to explain the mech- 



anism of the chirality dependence of the vortex pinning 
potential. We then present numerical result for the vor- 
tex pinning potential obtained from self-consistent order 
parameters. Our numerical result confirms the analyti- 
cal one. The chirality dependence of the vortex pinning 
would have influence on the hysteresis of the magneti- 
zation and the distribution of the magnetic field in sam- 
ples, which might be observed with SQUID and magneto- 
optical imaging techniques. 

To investigate the vortex pinning_we use the quasi- 
classical theory of superconductivity.EH We start with the 
Eilcnberger equation for the quasiclassical Green func- 
tion in the absence of the pinning, 



ffimt(«w n ,r,k) = -it: [ t _ 

? /imt Simt 



namely, 



iv-pk ■ Vpi m t + [iw n T z - A, g imt ] = 0, 



(1) 



(2) 



where the order parameter is A(r, k) = Uf x + 
if y )A(r, k) — (t x — if B )A*(r, k)l /2 and fi the Pauli ma- 
trices. The Eilenberger equation (p|) is supplemented 
by the normalization condition <?i m t(iw n , r, k) 2 = — 7r 2 l, 
and the commutator is [a, b] = db — ba. The vector 
r = (r cos <p, r sin </>) is the center of mass coordinate, and 
the unit vector k = (cos 9, sin 0) represents the relative 
coordinate of the Cooper pair. The cylindrical Fermi 
surface is assumed. We use maitSfin which K = feg = 1. 

Following Thuneberg et o/.jHdcil the effect of the pin- 
ning is introduced to the quasiclassical theory of super- 
conductivity as follows. The quasiclassical Green func- 
tion g in the presence of a point-like non-magnetic defect 
situated at r = R is obtained from the Eilenberger equa- 
tion 



iv F k ■ Vg + [iu n f z -A,g] = [i, g imt ]5(r') 
and the t matrix due to the defect 



(3) 



i(iu/ v 



v 
D 



i + N v(g imt (iuj n ,r' ,k)) g 



(4) 



= 1 



where r' = r R, the denominator D 
(irN v) 2 [(g im t)j + (/imt)e(/il nt }e], the average over the 
Fermi surface (• • ■)$ = J ■ ■ ■ dd/2%, the normal-state den- 
sity of states on the Fermi surface No, and we assume the 
s-wave scattering v when obtaining Eq. (^) . We define a 
parameter a = (7rA r w) 2 /[l + (7riVoi>) 2 ] , which measures 
how strong the scattering potential of the defect is. 
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The free energy in the-BEGsence of the defect is, at the 
temperature T, given asE3~E3 

5n(R) = N T^ dX^JdkJ drTi[Sg x A b ], (5) 

where 5g\ = g — <7; m t is evaluated at A = AAf,, and 
Aft is the order parameter in the absence of the defect. 
Equation (^) represents the difference in the free energy 
between the states with and without the defect, and then 
gives the vortex pinning potential <5f2(R). 

For the chiral p-wave pairing state d = z(k x + ik y ) = 
zexp(i#), it is known that the order parameter around a 
single vortex, A&(r, k) [= A&(r, <j>;9)], has two possible 
forms depending on whether the chimlito-and vorticity 
are parallel or antiparallel each other .B'tLZHia One form is 

A+- (r, <p; 9) = A + (r)e j(e ~^ + A_ (r)e l( - e+ ^ , (6) 

where the chirality and vorticity are antiparallel (Case 
I). The other is 

A+ + (r,<f>;6) = A+{r)e l{9+ ^ + A^{r)e l{ ~ d+ ^\ (7) 

where the chirality and vorticity are parallel (Case II). 
Here, the vortex center is situated at r = 0, the dominant 
component A + (r — > oc) = Abcs^), and the induced 
one A_(r — > oo) = 0. Because of axisymmetry of the 
system, we can take A±(r) to be real. 

First we analytically investigate the vortex pinning po- 
tential. We discuss the quantity 5tt(R = 0), where both 
the defect and the vortex center are situated just at the 
origin r = (R = |R|). From the quasiclassical view- 
point, the quasiparticles inside the vortex core, subject 
to the Andreev reflecti€«i|-*un along straight lines called 
as quasiparticle paths£3~E3 We consider the quasiparti- 
cle paths which go through the origin r = 0. On those 
paths with zero impact parameter, the position vector is 
parallel to the direction of the quasiparticle path (i.e., 
r || k), and therefore <fr — 6, 9 + tt. In this situation 
(<f> = 9), from Eqs. (JsJ) and (0), the order parameter on 
the path is 

A+-(r^ = 9;9) = A+(r) + A_(r) (8) 

in case I, and 

A++(r, 4> = 0; 9) = [A+(r) + A_(r)] e™ (9) 

in case II. The cancellation between the chirality and 
vorticity occurs in Eq. di) and not in Eq. (||). Of im- 
portance is the resultant difference in the phase factor of 
these order parameters. 

On the basis of an analysis of the so-called zero-core 
vortex model in Ref. ^3|, the matrix elements of— «; m t at 
the vortex center are approximately obtained asE3 

SWt = + lA] 2 ^ 1 , / imt = -Aw" 1 , // mt = A*w~ l , 

(10) 



where A = A+^r -> oo, = 9;9). Here, Eq. @ is 
obtained with assuming that the amplitude of the order 
parameter is constant (i.e., zero core) around the vortex, 
which is the only approximation in this analysis. Insert- 
ing the order parameter of Eq. (||) into Eq. (|l(]) , we ob- 
tain the anomalous Green functions integrated over the 
Fermi surface as, in case I, 

(/imt)« = /imt, (flat) 9 = /imt (H) 

because of the absence of any phase factors in Eq. (g), i.e., 
because of the cancellation between the chirality factor 
exp(i#) and the vorticity factor exp(— i<j>) in Eq. (^]). On 
the other hand, in case II, 

(firat)e = 0, (fL)0=0 (12) 

because of the phase factor exp(2i9) contained in Eq. (||). 
The diagonal component of (gi m t)e is (g\jat)e = Simt both 
in cases I and II. Consequently, in case L (g\ait)e — 9imt 
and we obtain [t,& m t] = from Eq. (Q). In case II, 
(ffimt)e |imt and [i,g imt ] ^ generally. 

When [i, g- im t] = 0, the Eilenberger equation (||) in the 
presence of the defect is identical to Eq. (||) (the equa- 
tion in the absence of the defect), namely, the defect has 
no influence on the Green function and the free energy. 
From this and the above results of the analysis of the 
factor [t, <7i m t]; we nn d that <5f2(0) = in case I when the 
chirality is antiparallel to the vorticity, and <5f2(0) ^ in 
case II when the sense of the chirality is the same as that 
of the vorticity. It means that the vortex pinning depends 
on the chirality in the chiral p-wave superconductor. 

The above analytical result is based on the zero-core 
vortex model, i.e., on the non-self-consistent (constant) 
amplitude of the order parameter. We next investigate 
the vortex pinning potential Sfl(R) numerically with the 
self-consistent order parameters around the vortex which 
have the forms of Eqs. (o) and (0). As self-consistent 
amplitude A±(r) in Eqs. ml) and (fa), we adopt numerical 
data which we have obtained in Ref. [l8] by solving self- 
consistently the Eilenberger equation. 

In Fig. [I], we show the numerical results for Sfl(R) in 
the Born limit [a 1) and the unitary limit (a — > 1). 
We present those results for the chiral p-wave pairing 
and the isotropic s-wave one. As noted in Fig. |l|, in the 
case of the s-wave pairing (dot-dashed lines), the dif- 
ference in the free energy between the states with and 
without the defect, 5Q(R), is equal to zero at R — > oo 
(R is the distance between the vortex center apd the de- 
fect). This is because the Anderson's theoremEJ is valid 
far away from the vortex core. On the other hand, in 
the chiral p-wave pairing cases (solid and dashed lines), 
<5fi(i? — > oo) is finite and positive as seen in Fig. [I]. The 
quantity Sfl(R — > oo) is equal to the loss of the conden- 
sation energy due to the pair breaking effect of the defect 
far away from the vortex core (i.e., the breakdown of the 
Anderson's theorem). As noted in Figs. |l|(a) and[l](b), at 
T = 0.8T C (high temperature), the condensation energy 
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loss in bulk Sfl(R — > oo) dominantly contributes to the 
depth of the vortex pinning potential 5Q(R), i.e., to the 
vortex pinning energy. From Figs. |l|(a) and |l](b) it is no- 
ticed that the vortex pinning energies of the chiral p-wave 
pairing cases at a high temperature are about 10 times 
larger than those of the s-wave pairing case. This en- 
hancement of the pinning effect is due to the breakdown 
of the Anderson's theorem, and then it must be .a .pom-, 
mon feature of unconventional superconductors. oEIrH 
For example, in the case of high-T c cuprates, this may 
be one of— the reasons why smalL-Boint defects such as 
Zn atomsE-3 and oxygen vacancicsE3 are efficient pinning 
centers. 

As noted in Figs. |(c) and|(d), at T = 0.2T C (low tem- 
perature), the contribution of the vortex core (R ~ 0) to 
the depth of 5il(R) is nonzero in case II (dashed lines). 
Here, the contribution of the vortex core means the en- 
ergy gain due to the presence of the scattering center in 
the vortex core. In contrast, the depth of 5fl(R), i.e., 
the vortex pinning energy, is determined in case I (solid 
lines) only by the loss of the condensation energy far 
away from the vortex core. It is noticeable that certainly 
5il(R — 0) equals to zero in case I. This numerical result 
confirms the analytical one discussed above. The vor- 
tex pinning energy depends on whether the chirality and 
vorticity are antiparallel (solid lines) or parallel (dashed 
lines). Especially in the Born limit, the difference in the 
vortex pinning energy is eminent as noticed in Fig. [|(c), 
because in this limit the loss of the condensation energy 
in bulk is relatively small compared to the contribution 
of the vortex core to the depth of Sfl(R). 

In general, the two chiral states of cases I and II can 
coexist as domain structure in samples under magnetic 
fields. The spatial gradient of the magnetic field in a 
sample is proportional to the local strength of the vor- 
tex pinning in the critical state. In terms of the present 
chirality-dependent vortex pinning, the gradient inside 
the domain of the case-II state is predicted to be steeper 
than that inside the domain of the case-I state. This may 
be experimentally observed as a signature of the chiral 
state. Also the domain structure of the two chiral states 
depends on the hysteresis of applied magnetic field, and 
therefore the present chirality-dependent vortex pinning 
may affect the hysteresis curve, during multiple cycles of 
the magnetization as observedEl in UPt3. 

In the case of the usual winding-1 vortex A oc exp(i0), 
the chiral "p-wave" pairing k x ± ik y = exp(±i#) is essen- 
tial for the cancellation between the chirality and vortic- 
ity. If winding-2 vortices A cx exp(2i<^>) are realized in 
a chiral <i-wave state k x — k y ± ik x k y = exp(±2i9), the 
same kind of cancellation occurs. 

We comment on the relation of the present vortex pin- 
ning to the superconducting gap structure in Sr 2 Ru04. 
In this material, it has been pointed out from experi- 
ments that-the gap has line nodesH~H and little in-plane 
anisotropy.Q Models for the gap structure consistent with 
those experimental facts were proposed, in which there 
existed horizontal line nodes perpendicular to the axis 



of the cylindrical Fermi surface.™} Now, for the present 
theory of the chirality-dependent vortex pinning, what 
is important is that the Fermi surface averages of the 
anomalous Green functions (i.e., the average of the or- 
der parameter except for the chiral part) are finite as in 
Eq. ( pr|) . The present chirality dependence of the vortex 
pinning does not occur, if the order parameters have sign 
changes on all Fermi surfaces relevant to superconduc- 
tivity as A(k) ~ exp(±i0) cos(cfc z ). It occurs, if there 
are no sign changes as A(k) ~ exp(±i0)\ cos(cfc z )|. In 
another case,B the chirality dependence is expected to 
occur when the order parameter is nodeless on the ma- 
jor Fermi surface with dominant density of states, even if 
there are gap nodes and sign changes on the other minor 
Fermi surfaces. 

In conclusion, we investigated the elementary vortex 
pinning potential dfl(R) on the basis of the quasiclas- 
sical theory of superconductivity. In the chiral p-wave 
pairing state, Stt(R) was dependent on the sense of the 
chirality relative to the vorticity at a low temperature. In 
terms of the present chirality-dependent vortex pinning, 
a theoretical analysis for anomalies in the hysteresis pf 
the magnetization observed experimentally in S^RuCUEa 
would be interesting and is left for future work. 
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FIG. 1. The vortex pinning potential as a function of the 
distance R between the vortex center and the defect. Solid 
lines [p(H — )] correspond to the case of the p-wave pair- 
ing with the chirality antiparallel to the vorticity (Case I). 
Dashed lines [p(++)] correspond to the case of the p-wave 
pairing with the chirality parallel to the vorticity (Case II). 
Dot-dashed lines correspond to the case of the isotropic 
s-wave pairing. T c is the superconducting critical tempera- 
ture. The distance R is normalized with the coherence length 
£o = «f/A B cs(T = 0). 
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